Characterizing dynamic length scales in glass-forming liquids 
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In their correspondence [1] Flenner and Szamel (FS) compare the temperature dependence of an alternative dynamic 
length scale, £4, with the one of ^ dyn that we have studied in [2]. Using computer simulations of the same system, they 
conclude that these two length scales have a different temperature dependence. In particular, the former does not follow 
the striking nonmonotonic temperature dependence we reported for the latter. While both types of measurements 
aim at quantifying the spatial extent of dynamic correlations in supercooled liquids, the two procedures differ on 
essential points, which we now discuss. 
I . While we measured up to what distance the value of the relaxation time is affected by the presence of an amorphous 
^C) ' wall [2], FS quantify instead the spatial extent of spontaneous dynamic fluctuations at low wave- vectors for a fixed 
timescale (the bulk relaxation time). These two measurements need not be directly related, although they seem 
to coincide at moderate temperatures (see Fig. la of FS). While four-point functions as measured by FS have 
played a pivotal role in previous analysis of dynamic heterogeneity, theoretical work has also revealed a number of 
shortcomings. Most notably, four-point functions display a strong dependence on the statistical ensemble chosen to 
perform measurements. As a result, they receive distinct contributions from density and energy fluctuations and 
show, even in idealized cases, complex scaling properties [3] which complicates the direct extraction of a correlation 
length. A second difficulty lies in the fact that these various contributions have different temperature dependences, 
with a crossover taking place very close to the mode-coupling temperature where the subtle effects reported in [3] 
occur. We remark that this corresponds to the temperature scale where £4 and t; dvn start to differ. These known 
weaknesses of four-point functions had in fact motivated our study of an alternative correlation length scale that is 
free of such ambiguities. 

Given these important differences, it is not clear how a nonmonotonic temperature evolution of dynamic correlations 
will manifest itself in the numerical data of FS. While FS argue that the evolution of £4 with the relaxation time 
l— ~" presents a crossover (Fig.lc in FS), we point out that the shown data does not clearly show such a change in behavior. 
Another possibility, not explored by FS, could be that the functional form of the correlator S^(q,t) changes with 
■ temperature, in agreement with the idea that the geometry of the relaxing entities changes with temperature, as 
t^J- I argued in [2]. Since the reported effect is small (Fig. 2b in [2]), one would presumably need a relative accuracy of 
t-H . Si(q,t) of better than 1% at low wave-vectors. The data shown by FS demonstrates that at present this remains 
technically extremely difficult. 

. Finally, the crossover we report also coincides with the emergence of nontrivial static correlations in the system, a 
result that cannot be obtained using purely dynamic correlations. Overall, this suggests that the approach of pinning 
, particles, for which no a priori knowledge of the nature of the relaxing entities is needed, is extremely helpful to 
study small but important effects in the relaxation properties of glass-forming systems that cannot be drawn from 
the measurements of FS. We point out that since the publication of [2], we independently confirmed the qualitative 
change of transport properties in our system near the mode-coupling crossover by analyzing the temperature evolution 
of dynamical finite size effects [4] . 
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